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Abstract 

■ The dependence of the fractal behaviors of the pomeron induced system in deep 

■ inelastic lepton-nucleon scattering upon the diffractive kinematic variables is found 

■ rather robust and not sensitive to the distinct parameterization of the pomeron flux 

■ factor and structure function. A feasible experimental test of the phenomenological 

o ■ 

00 ' pomeron-exchanged model based on the fractal measurement in DESY ep collider 

r~| ■ HERA is proposed. 
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Diffractive processes in high energy colhsions have long been described by using the phe- 
nomenology of Regge theoryS by the t-channel exchange of mesons and, at high energy, by 
the leading vacuum singularity, i.e. the pomeroni. Many calculations based on the pomeron 
theory has been pursued concerning with the collective aspects of the diffractive processes, 
such as cross section of hard diffractioni^, the distribution of large rapidity gapl^B, jet 
productionii'i in hard diffractive processes etc.. But because of the ignorance for the na- 
ture of the pomeron and its reaction mechanisms, there exists different kinds of approaches 
and parameterizations of pomeron in current Regge theory. The calculation results in the 
above-mentioned literature were found very sensitive to the distinct parameterization of the 
pomeron. So the experimental data concerning those features can constantly be conformed 
by pomeron model with tunning all these unknown parameterizations, although the param- 
eterizations are quite different for different respects of dataf], which makes it rather difficult 
to conclude whether pomeron theory indeed works in the diffractive process. In this respect, 
it is nature to ask the following questions: Is there a way to test and justify the pomeron 
exchange model by using current diffractive experiment equipment (e.g. DESY ep collider 
HERA) while the criterion to the measurements do not depend upon concrete parameter- 
ization of pomeron? If yes, what is the characteristic behavior of the pomeron exchanged 
model in the expected experimental measurements? 

Having in mind that the fractal and fluctuation pattern of the multiparticle production 
reveals the nature of the correlations of the spatial-temper evolutions in both levels of 



^For example, the earlier data of UA8 ColIaborationEl about jet production in CERN SppS- 
Collider dictated a soft gluon content of pomeron, while newer and more comprehensive data of 
UA8 Collaborationii required a hard pomeron content; the analysis of structure function by ZEUS 
Collaboration0 in DESY ep collider HERA needed both hard and soft component of the pomeron; 
the recent study of structure function of QCD evolution by HI collaboration0 favored a rather 
peculiar one-hard-gluon pomeron. 



1 



parton and hadron and is, therefore, sensitive to the interact dynamics of the high-energy 
processEl'El, it has been proposed to investigate the fractal behavior of the diffractively 
produced system by calculating the scaled factorial moments of the multihadronic final 
stated. In this note, we find the dependence of the fractal behavior of the pomeron induced 
system upon the diffractive kinematic variables is rather robust and not sensitive to the 
different parameterization of the phenomenology pomeron model in the deep inelastic lepton- 
nucleon scattering (DIS). So the characteristic plot about fractal behavior of the diffractively 
produced system can be considerated as a clear experimental test of the pomeron exchanged 
model by DESY ep collider HERA. 

The fractal (or intermittency) behavior of the diffractively produced system in DIS (and 
also hadron-hadron collider) can be extracted by measuring the g-order scaled factorial 
moments (FMs) of the final-state hadrons excluding the intact proton from the incident 
beam, which are defined bylll 

p / c N 1 ^ {nmirim - 1) . . . {n„, -q + l)) . . 

where, x is some phase space variable of the multihadronic final-state, e.g. (pseudo-) rapidity, 
the scale of phase space Sx = Ax/M is the bin width for a M-partition of the region Ax 
in consideration, Um is the multiplicity of diffractively produced hadrons in the mth bin, 
(■ ■ ■) denotes the vertical averaging with the different events for a fixed scale 6x. Since the 
factorial moments Fg can rule out the statistic noise around probability Pm for a particle to 
be produced in the small phase space of final state, and associated directly with the scaled 
probability moments 

1 (v'i ) 

it is clear that Fg would saturate to some constant with decreasing 6x to the some typical size 
6xo (e.g. the correlations length from resonance decays) if there were no dynamic fluctuation 
in the multiparticle producing process and the probability distribution of the phase-space 
densities were smooth, especially in small scale of the phase space. The manifestation of the 



fractality and intermittency in high energy multiparticle production refers to the anomalous 
scahng behavior of FmIiIO 

Fg{5x) ~ {Sxy^" ^ as M ^ cx), 5x 0. (3) 

The g-order intermittency index (pq can be connected with the anomalous fractal dimension 
dq of rank q of spatial-temporal evolution of high energy collisions asEl 

dq = (t>q/{q - 1) (4) 

Pomeron exchange in Regge theory has be used to describe successfully the main features 
of the high energy elastic and diffractive processlll. While waiting for the experimental 
measurements of DESY ep collider HERA (and also hadron-hadron collider) for the above- 
mentioned fractal behavior of the diffractively producted system, let us now take a closer 
look at what we can learn from the current pomeron exchange theory. 

The deep inelastic diffractive lepton-nucleon scattering has been well modeled by several 
Monte Carlo generators such as POMPYtIiiI and RAPGApi, which are based on the 
sumption that the process can be considered as taking place in two stepsi. Firstly, a b 
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hadron emits a pomeron (P) with longitudinal momentum fraction xp; in the second step 

this pomeron interacts with the other beam lepton in a large momentum transfer process 

between the basic constituents of the pomeron and the lepton. This pomeron factorization 

allows the diffractive hard scattering cross-section to be written as 

d^a{ep ^e + p + X) _ , . d'^crjeV ^ e + X) 

dxpdtdpdQ^ -/pPl*,a;pj ^^^^2 ' 

where —Q^ and t are the mass-squared of virtual photon and exchanged pomeron respec- 
tively, and (3 = Xb/xp the momentum fraction of parton in pomeron. The first factor in 
Eq. is the pomeron flux, i.e. probability of emitting a pomeron from the proton. For 
lack of the understanding for the reaction mechanism of the pomeron, there exist different 
approaches for the parameterization of the pomeron flux: 

In standard Regge theory with a supercritical pomeron trajectory, ap(t) = 1 + e + a'pt, 
the pomeron flux was given asllli 



f \ /^pp(^) l-2ap(t) /„\ 

iOTT 

where the couphng of the pomeron to the proton, /?pp(t) was parametrizedi as /5pp(t) = 
6.3e^'^^* according to the t-dependence of cross section in higher xp-range of ISR datali. 
By approximating the cross section of proton and pomeron with a constant (Xpp = Imb, 
Ingelman and Schleini compared the factorization expression of hard diffractive pp process 
with the single diffractive differential cross section at the SPS collider^. They obtained the 
pomeron flux in a proton as 

/pp(t,xp) = TT^ — (6.38e-«l*l + 0.4246-=^!*'). (7) 

Z.o Xp 

Alternatively, Donnachie and Landshoff assumed the pomeron-photon analogy in the couples 
to quark in nucleon and gave the pomeron flux a 



fMt,x^) = ^-^m)]%-"'''^'^ (8) 

with coupling constant of the pomeron to a quark /J^p = 3.24GeV~^ and the elastic form 
factor 

F (f\ _ 4mg-2.79f 1 / X 

The second factor in Eq. (^, i.e. lepton-pomeron hard cross section, can be calculated 
in the way that 

rfV(eP^e + X)_ /■ (iVhard(e + parton->e + X) 

if we could figure out the density G{f3) of the quarks and gluons with fraction (3 of the 
pomeron momentum, needless to say, hard scattering cross section (icrhard(e+parton e+X) 
is to be computed according to standard QED and the perturbative QCD correction. But 
the pomeron structure function is still a main uncertainty in pomeron model and it is 
even unknown whether the pomeron consist mainly of gluonsil'i'iii or of quarklil, although 
measurements of hard diffractive scattering have been performed in both lepton-hadron and 



hadron-hadron collider (see Footnote '''). Two extreme gluon densities has been extensively 
used in Ref.l'§i'i0 and many other literature, i.e. 

(3G{P) = 6/3(1-/3), (11) 



PG{(5) = 6(1 - (3)\ (12) 

The first one corresponds to the "unrealistically hard" gluon distribution, where two gluons 
share the pomeron momentum; in latter case, the gluons in the pomeron are as soft as that 
in the proton. Taking into account the information from Regge theory in /3 — > and the 
power counting rules in f3 1, Berger et al.i parametrized the gluon distribution in the 
pomeron as 

(3G{/3) = (0.18 + 5.46/3)(l - /?). (13) 

In above three parametrizations of structure function of pomeron, the normalization of all 
parton distributions are, by default, chosen to fulfill the momentum sum rule 

^ dppCiP) = 1. (14) 

But it is not clear that this relation must hold for the pomeron which is a virtual exchanged 
object that need not behave as a normal hadron state. In the approach of Donnachie and 
Landshoffi, the dominating quark density was 

PG{(3) = -C7cf3{l - (3), with C = 0.23, (15) 
3 

the flavors summation and the momentum integral of which is a factor 7.5 lower than the 
normalization in Eq. (p!4D. 

In addition to these theoretical uncertainties there is also a uncertainty in the evo- 
lution of the parton densities of the pomeron. Numerical calculations using ordinary QCD 
evolution equation (Altarelli-Parisi or DGLApi), and GLR-MQi equation in which the 
inverse recombination processes of partons has been taken into account, turned out thatJ^l 



the evolution of the pomeron structure function can be very much different depending 
upon whether the non-hnear recombination term of the QCD evolution equation is included 
or not. Furthermore, depending upon the initial parton distribution at a given momentum 
scale which is unknown, the size of nonlinear term may become too large for the QCD evo- 
lution equation to be reliable without further, but also unknown, correction. By assuming 
both leading and subleading Regge trajectory with a flux akin to Eq. ( ||), a fit according 
to the NLO DGLAP evolution equations to HERA data0 of F.^^^\xp, P,Q'^) has favored a 
rather peculiar "one-hard-gluon" distribution for the pomeron (see also Footnote ^). Since 
what we try to pursue in this note is to find out whether and in which range the fractal 
behavior of pomeron induced system depend on varieties of different parameterization of 
pomeron, we leave the possible anomalous scaling behavior in the QCD evolution processes 
to the further discussionll. 

In typical kinematic region of hard diffractive processes of DESY ep collider HERA (say, 
i.e. Mx > l.lGeV, and xp < 0.1), majority properties of the diffractive events can be 
well reproduced by RAPGAP generator (see, e.g. In the following intermittency 

analysis of the lepton-nucleon diffractive process, we use RAPGAP generator^ to simulate 
the pomeron exchange processes, in which the virtual photon (7*) will interact directly 
with a parton constituent of the pomeron for a chosen pomeron flux and structure func- 
tion. In addition to the 0{aem) quark-parton model diagram (7*5 q), the photon-gluon 
fusion (7*(7 qq) and QCD-Compton (7*5 qg) processes are generated according to 
the 0{aemOts) niatrix elements. Higher order QCD corrections are provided by the colour 
dipole model as implemented in (ARIADNE)il, and the hadronization is performed using 
the JETSET0. The QED radiative processes are included via an interface to the program 
HERACLEsi. 

By chosing the pomeron flux factor /pp(t, Xp) and the pomeron structure function G{I3) 
as shown in Eq. (|]) and Eq. ( pA]) respectively, we generate 100,000 MC events, and cal- 
culate the second-order factorial moments in 3-dimensional (?7,p^,0) phase space, where 
the pseudorapidity 77, transverse momentum p± and the azimuthal angle are defined with 



respect to the sphericity axis of the event. The cumulative variables X translated from 



were used to rule out the enhancement of FMs from a non-uniform inclusive spectrum 



decreasing scale of the phase space is shown in Fig. 1(a) in double logarithm. There exists 
obviously anomalous scaling behavior in the pomeron induced interaction, so we fit the 
points in Fig. 1(a) to Eq. (^ with least square method and obtain the intermittency index 
02 = 0.428 ±0.003. In Fig. 1, we have also shown the Monte Carlo result of the second-order 
factorial moments for different kinds of the parameterization of the pomeron flux factor and 
structure function. In Fig. 1(a), (b), and (c), we keep the pomeron structure function G{P) 
fixed but vary the pomeron flux factor /pp(t,xp) as Eq. (^, (^, and (H) respectively. The 
fractal behaviors of the pomeron induced system keep almost unchanged for different flux 
factor. On the contrary we keep the pomeron flux factor fixed in Fig. 1(a), (d), (e), and (f), 
but vary the pomeron structure function as Eq. (pA]), (p!^), (|13]), and (|T^) respectively. For a 
given pomeron flux, the fractal behaviors become weaker when the pomeron become softer. 
In Fig. 1(d) the parton distribution is as soft as that in proton, the intermittency index is 
smallest, which is understandable since if the hard parton in pomeron is involved it is more 
possible to evoke jets and then the anomalous short-range correlation in the final-state so 
that the intermittency index increases, and vice versa. 

It is of the special interest to investigate the dependence of the fractal behavior of the 
pomeron induced system upon the diffractive kinematic variables. We generate 500,000 
events by RAPGAP Monte Carlo generator, and divided the whole sample into 10 sub- 
samples according to the diffractive kinematic variables, e.g. xb- For each subsample, 
we calculate the second order scaled FM and the intermittency index 02, and to see how 
the fractal behavior of the pomeron induced multihadronic final states depends upon the 
considerated kinematic variable. In Fig. 2 is shown the dependence of the second order in- 




(16) 



p{x) of the final produced particleg22l. The obtained result of second-order FM 




versus the 
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termittency index 02 on the different diffractive kinematic variables. Since it is well known 
that the gluon density increase sharply as xb decreases in small-x^ regionlll, the MC result 
from pomeron model in Fig. 2(a) means that the anomalous fractal dimension d2 of the 
diffractively producted system decrease with increasing gluon density, which is not incon- 
ceivable if one take into account the fact (see, e.g.liiEi) that the effect of superposition of 
fractal systems can remarkably weaken the intermittency of whole system. Obvious depen- 
dence of 02 on pomeron momentum fraction xp of a hadron and parton momentum fraction 
/? of a pomeron as shown in Fig. 2(b) and (c) implies that, the intermittency calculated 
here can not be only referred to the hadronization processes and there should be substan- 
tial correlations between the fractal calculation and pomeron dynamics. In Fig. 2(a) and 
(b), the intermittency index 02 is less than for the lower xb and xp, which can be im- 
puted to the constraint of the momentum conservation in the high energy procesd^. Since 
the Leading Proton Spectrometer (LPS) has been used in ZEUS detector to detect protons 
scattered at very small angles (say, < 1 mrad), which make it possible to measure precisely 
the square of the four-momentum transfer t at the proton vertex, we also showed in Fig. 
2(e) the t-dependence of second order intermittency index in the t-region of LPS detector, 
i.e. 0.07 < —t < 0.4GeV^. To be different from the results of other kinematic variables, the 
fractal index for the 7*P system doesn't depend upon the t. 

Especially, we calculate the intermittency index shown in Fig. 2 using different kinds 
of pomeron parameterization. We denote the different shapes of points in Fig. 2 for the 
different kinds of the pomeron parameterization, just in the same way as that in Fig. 1. 
In conventional investigation, the pomeron theory has been used to compare with the data 
about cross section of hard diffract ionlHi, the rapidity distribution of large rapidity gapifl, 
and jet rapidity distributionii and jet shapeBi etc., where the results of the pomeron model 
in such quantities concerning with the collective nature of diffractive process were found 
very sensitive to the parameterization of the pomeron, and the experimental data in dif- 
ferent aspects preferred different kinds of parameterizationi^i3 (see also Footnote It 
is remarkable that the dependence of the intermittency index, which concerned with the 
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inherent scaling behaviors of diffractive processes, upon the diffractive kinematic variables 
in this implementation of the pomeron exchanged model are rather robust and almost the 
same for the different parameterization of the pomeron flux and structure function! 

In conclusion, we have given an outline of the mainly uncertainty of the pomeron ex- 
changed model on the distribution functions for finding partons in a pomeron and for finding 
pomeron in a hadron, which is confronted in the conventional hard diffractive calculation 
(see, e.g. Re^^0). By the intermittency analysis according to Monte Carlo implementa- 
tion of RAPGAP we present a characteristic plot for the pomeron exchange model, which is 
independent of aboved- mentioned uncertainty. So in order to test and justify the pomeron 
theory to the diffractive processes, it is urgent and feasible to check this characteristic frac- 
tal plot in DESY ep collider HERA. And the substantial revision would be necessary in 
the manner in which we have treated diffraction by using the current pomeron theory if it 
should turn out that experimental measurements differ drastically from this characteristic 
plot presented here. 
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FIGURES 

Fig. 1. The second-order scaled factorial moments F2 getting from MC simulation of RAP GAP 
generatori versus the number M of subintervals of 3-dimensional {'i],p±,4>) phase space in log-log 
plot, and the intermittency index (/>2 correspondingly. The different kinds of points denote different 
parameterization of the pomeron flux factor fpp{t.xp) and structure function G(/3), i.e. /pp(t.a;p) 
keeps fixed as Eq. ( pT| ) but G{(3) varies as Eq. (|6|), (|7|), and (^) in Figs, (a), (b), and (c) respectively; 
on the contrary, G{(3) keeps fixed as Eq. (^) but /pp(t.xp) varies as Eq. (|ll|), (|l2|), (13), and (15) 
in Figs, (a), (d), (e), and (f) respectively. 



Fig. 2. The dependence of second-order intermittency index (j)2 in RAPGAP0 Monte Carlo 
implementation upon different kinematic variables. The different shapes of points denote different 
parameterization of pomeron flux factors and the structure functions in the same way as that in 
Fig. 1. 
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